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Abstract 

We propose a fundamental relation for a classical ideal gas that is valid at all temper- 
atures with remarkable accuracy. All thermodynamical properties of classical ideal 
gases can be deduced from this relation at arbitrary temperature. 

The famous equation of state for an ideal classical gas is 

PV = Nk B T. (1) 

Interestingly, a classical gas obeys this relation at all temperatures as long as it is ideal, 
i.e., the Hamiltonian of the system does not depend on the co-ordinates of the particles at 
all. The proof is simple, and appears in most textbooks We will also provide the proof 
later. 

An equation of state, however, does not specify a system completely 0. For example, 
from the equation of state in Eq. (|l|), we cannot find the entropy of the system, and 
many other properties for that matter. Of course if we have a fundamental relation for the 
system, it contains all thermodynamic information about the system including the equations 
of state themselves [0. These are, for example, relations of the type S = S(U,V,N) or 
U = U(S,V, N), which express the entropy or the internal energy as a function of other 
extensive parameters of the system. Legendre transforms of these equations work just 
as well, like the Helmholtz free energy A as a function of T, V and N. However, it 
is often difficult to obtain such relations in closed forms which would be valid for any 
temperature. The Sackur- Tetrode relation, for example, is a fundamental relation of the 
form S = S(U, V, N), but unlike the relation in Eq. (|l|), it is valid only if the gas is non- 
relativistic, i.e., if the temperature is small in the sense that (3mc 2 3> 1, where m is the 
mass of the gas particles. Our aim in this article is to suggest a fundamental relation for 
the classical ideal gas that can be used at any temperature. 

Since the gas is assumed to be ideal, the energy of any particle in the gas depends only 
on its momentum. At a momentum p, let us denote the energy of a particle by s(p). The 
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single-particle partition function in the canonical ensemble is then given by 

where g denotes the degeneracy due to internal degrees of freedom, A is the Compton 
wavelength of the particle divided by 2tt, 

X = —, (3) 

mc 

and /(/3) is a dimensionless integral defined by 

/0S) = 7 — ^ / dpfe-M. (4) 
[mcy Jo 

The canonical partition function for a system of N particles is then given by 

7 _Q N 

with the Gibbs correction factor for identical particles. This gives 



(5) 



\nZ = N 



In V + In / - In iV + 1 + \n(g /2% 2 X 3 ) , (6) 



where we have used Stirling's formula for IniV!. It is to be noted that f(/3) is the only 
component which is not exactly known at this point. Our aim would be to determine it. 

For this, we start by evaluating some thermodynamic quantities. The pressure of the 
gas is given by 

l/£taZ\ = N_ 

irrespective of the functional form of f{/3). This is Eq. (ffl), and this is why it is valid for 
any temperature. 

Next we look at the energy density. This is given by 

1 (d\nZ\ Ndlnf 
6 V{ d(5 ) VN V dp ■ W 

Thus, if we define a quantity w by 

P = WQ, (9) 

Eqs. (0) and (|) show that 

w ^ dp d In (3 ^ 
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Figure 1: The ratio Pj g for different temperatures. The solid line is obtained by numerical 
integrations from Eq. (|13|), and the points correspond to the fit of Eq. (|l^). The dashed 
horizontal lines are drawn at |, | and |. 



Using the definition of f(/3) from Eq. (|4]), we can then write 



dp p 2 e- pe{p) 

w = ^ -. (11) 

dpp 2 pe{p)e-P £{p) 



It is easy to see that if e{p) oc p x , this formula gives w = A/3, which gives the correct results 
for the very low and very high temperature limits [|TJ. 

For arbitrary temperatures, we should use the relativistically correct formula for e(p). 
For an ideal gas, there is no potential energy. The kinetic energy of a particle with momen- 
tum p is given by 



e(p) = \Jm?c 4 + p 2 c 2 — mc 2 . (12) 
Putting this into Eq. ([□]) and making a change of variable, we can write 
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W= ,oo J ° - ; . (13) 



where ^ is a dimensionless variable for the inverse temperature: 

£ = [3mc 2 . (14) 

The integrations can be performed numerically for any value of £. The results are shown 
by the solid line in Fig. [l|. 

The shape of the curve is very similar to a hyperbolic tangent curve. We can try a fit 
of the form 
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i + itanh(aln(e/6)), (15) 



with two parameters a and b which are both positive. This has the correct behavior for 
both extremes: w = 1/3 for very high temperatures (log£ — > — oo) and w = 2/3 for very 
low temperatures (log£ — > +00). For reasons described below, we must take a — |, which 
means that our fit is of the form 

«, = I + -|^. (16) 
2 6 £ + b v ; 

The choice for the parameter b will be made later, in Eq. (p9|). However, in order to give 
a preview of the goodness of this fit, we anticipate the value and show the resulting fit as 
points in Fig. [I]. As the figure suggests, the fit is very good for all temperatures. In fact, 
the agreement between the fit and direct numerical calculation is considerably better than 
1% everywhere. 

It has to be remarked that e(p), and consequently g, do not include the contribution 
from the rest mass energy of the particles. However, it is trivial to obtain a relation between 
P and the energy density p which includes the mass energy as well. Obviously, 

p — g + mc 2 n , (17) 

where n = N/V is the number density of particles. Using Eq. ([]]), we can write it as 

p = Q + tP (18) 

for a classical ideal gas. If we then define an equation of state in the form 

P = up, (19) 

the constant u will have a very simple relation with w which can be obtained by dividing 
both sides of Eq. (0) by P: 

- = - + £• (20) 

to w 



We now go back to the fit of Eq. fll6|), and notice that from Eq. (|T0"D, we can now 
determine In/. In fact, using the fit, we can write 

noting that din (3 = cHn£, since mc 2 is constant. Thus, 

m/ = -3[ d 4 1±IL 

= K-31n£ + ^hi(2£ + &), (22) 

where K is an integration constant. 

To determine this constant, we may use of the form for In / for large £. In this case, 
small momentum values dominate the integrand, so that we can approximate e(p) by p 2 /2m. 
Once this is done, the integration in Eq. (|J) can be exactly performed and one obtains 

1 t-3/2 



/(o^^r- (23) 
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On the other hand, for £ ^> 1, Eq. ( p2| ) gives 

ln/(£>>l)«#-^ln£+^ln2. (24) 

Notice that this has the correct £-dependece as the non-relativistic formula of Eq. ( p3f) . 
This is the reason why we had to take a = | in Eq. 015]). Moreover, this also allows us to 
determine .fT, viz., 

K = ilnvr - 21n2. (25) 

Putting this value of K into Eq. (|2~2"D, we obtain 

1 3 
In/ = -In7r-21n2-31nf + -ln(2f + 6). (26) 

The remaining constant h can be determined by using the asymptotic value of / for 
small £. In this case, large momentum values dominate the integrand so that we can put 
e(p) = cp. Once this approximation is made, the integration in Eq. (^) can be performed 
analytically and we obtain 

/(£«l)^2r 3 . (27) 

On the other hand, Eq. (^) gives 

ln/(£< 1) « -Invr-21n2-31n£ + -ln6. (28) 
2 2 

Comparing the two expressions, we obtain 

' I < 29) 

This finally completes the fit that we sought for. The fit is shown in Fig. 0, along with the 
values obtained by direct numerical integration. 

We would like to make a comment about our choice of the parameter b. We have chosen 
it so as to fit the low-£ end perfectly. However, the resulting fit is worst (though still better 
than 1%) for intermediate values of £, which shows clearly in Fig. [I]. If we take a slightly 
smaller value of £, the middle part fits much better. Of course it also means that the 
agreement for small £ gets a bit worse. To be precise, if instead of the value b given in Eq. 
(^9|) we choose some other value b' for the parameter under consideration, the values of In / 
mismatches by an amount |ln(&/6') for ( < 1. But since the value of /(/?) itself is very 
large for small values of £, the fractional error is very small. A choice like b = e, the base 
of natural logarithms, makes the overall fit much better, to better than 0.5%, for the entire 
range of temperatures. If we go further down, say b = 8/3, the fit becomes very good in 
the middle but considerably worse at the ends. 

We can put our fit for /(£) into Eq. (R|). This gives 



\nZ = N 




(30) 
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Figure 2: The integral f(/3) defined in Eq. @, plotted as a function of £ = f3mc 2 . The two 
dashed lines correspond to the asymptotic forms for f(/3) given in Eqs. fl23|) and (p7|) . The 
points are our fit from Eq. fl26|) , with the value of 6 taken from Eq. (p9|). 



where 6 is given in Eq. (|29|). The Helmholtz free energy of the classical ideal gas is then 
given by 




(31) 



If we recall the definition of £ from Eq. (|T4T), we realize that this is an expression for the 
free energy as a function of the volume, temperature and total number of particles. This 
is thus a fundamental relation which can be used at any temperature. All thermodynamic 
properties of an ideal gas can then be deduced from this relation at arbitrary temperature. 

This is the result. Of course it is not exact; it is a fit to the numerical results. But it 
is a very good fit, and the usefulness of such an expression cannot be overemphasized. In 
most practical situations, a closed analytic expression is much easier to use than numerical 
integrations. Our expression for the Helmholtz free energy in Eq. (|3lD provides such an 
analytic expression which can be used at any temperature to obtain all thermodynamic 
properties of a classical ideal gas with suitable manipulations. 

We can of course write the fundamental relation in alternative but equivalent ways. For 
example, we can write the fundamental relation for the entropy S by expressing it as a 
function of U, V and N. For this, we use Eq. fl3l|) to obtain 



S 



Nkf 




V,N 




(32) 



At this point, the right hand side is a function of V, N and T (through £). To put it in 
the form of a fundamental relation, we must eliminate £ and bring in the internal energy 
U. This can be done through the fit of Eq. (|30|) , which gives 

Recalling the definition of £ from Eq. ([141) , it can be easily checked that this gives U = 
3N/2/3 in the non-relativistic limit (£ ^> 1), and U = ?>N/(3 in the ultra-relativistic limit 
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(£ <C 1). In the general case, we can solve Eq. ( |33"D to obtain 



3Nmc 2 -bU + ^(3Nmc 2 - bU) 2 + 2AbUNmc 2 

4T7 



(34) 



Putting this back into Eq. fl32|), we obtain the expression for entropy as a function of the 
external variables U, V and N. As stated before, this is also a fundamental relation. The 
non-relativistic limit should now be acknowledged as U <C Nmc 2 , where Eq. ( |3"4"D gives 
£ = 3Nmc 2 /2U, which is large. Putting this solution into Eq. fl32|), we obtain 



S 



Nk 



B 



Ia^A 3 ; 2 WirNmc 1 ) 2 



In 



m 



3/2 



AT5/2 



+ 



2 ' 



(35) 



which is the Sackur- Tetrode relation 0. In the other extreme, when U ^> Nmc 2 and so 
£ = 3Nmc 2 /U < 1, we obtain 
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A^ 



In 



7T ^\/f/ 3 



(3tt^c) 3 A^ 4 



+ 4 



(36) 



if we use the value of b from Eq. (|29|) . This result can also be obtained analytically by using 
e = cp. 

I thank R. K. Moitra for the suggestion of including the expression for entropy in the 
paper. 
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